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SUMMARY 

This  report  summarizes  my  research  on  spectral  characteristics  of  excitation  and 
propagation  of  Pn  and  Lg  in  the  past  year.  Pn  and  Lg  spectra  from  the  last  eight  Lop  Nor 
explosions  and  many  nearby  earthquakes  are  collected  from  many  broadband  IRIS, 
CDSN,  Kyrghizstan  and  Kazakhstan  network  stations,  and  analyzed  to  obtain  source 
spectral  characteristics  and  path  Q.  Across  the  ten-station  Kyrghizstan  network,  time 
domain  Pn  amplitude  varies  by  as  much  as  a  factor  of  19,  despite  the  small  aperture  (< 
200  km)  of  the  network.  Spectral  ratios  of  large,  with  respect  to  small,  explosions  are 
associated  with  a  spectral  overshoot  that  is  more  significant  than  predicted  by  existing 
theoretical  models.  Pn  from  earthquakes  is  further  complicated  by  effects  of  non-isotropic 
radiation  patterns.  Several  efforts  are  made  to  reduce  the  effect  of  these  complexities  in 
Pn  spectral  analysis.  These  include  (a)  the  use  of  a  moving-window  average  of  spectra  of 
Pn  and  its  early  coda,  (b)  a  modification  to  the  inverse  method  of  Xie  (1993)  to  tackle 
problems  in  Pn  spectral  inversion;  (c)  stacking  spectra  over  similar,  large  explosions  to 
obtain  robust  estimates  of  path  Q;  and  (d)  the  use  of  multiple  constraints  on  source  spec¬ 
tra  and  path  attenuation  by  combining  Pn  spectral  inversion  with  analysis  of  inter-event 
and  inter-station  spectral  ratios.  The  main  findings  of  the  analysis  are:  (1)  for  both  explo¬ 
sions  and  earthquakes,  the  logarithm  of  seismic  moments  (Mq),  estimated  by  inverting  Pn 
and  Lg  spectra,  scales  with  body  wave  magnitude  (m^,)  linearly,  with  slope  of  about  1.0, 
(2)  for  earthquakes,  the  Mq  estimated  using  Pn  (Pn  Mq)  is  similar  to  Lg  Mq  at  a  given 
level,  whereas  for  explosions,  Pn  Mq  tends  to  be  higher  than  Lg  Mq;  (3)  for  both  Pn  and 
Lg,  Mo  scales  with  /“®,  with  a  being  close  to  4  for  explosions,  and  3  for  earthquakes;  (4) 
for  both  explosions  and  earthquakes,  at  the  same  Mq  level,  Pn  /<-  are  much  higher  than  Lg 
/c  (by  a  factor  of  about  4  for  earthquakes  and  5  for  earthquakes,  respectively),  (5)  for 
explosions,  the  significant  spectral  overshoot  in  Pn  causes  the  Pn/Lg  ratio  to  reach  maxi¬ 
mum  around  Pn  f^.  This  maximum  is  the  dominant  cause  for  that  ratio  to  be  significantly 
higher  for  explosions  than  for  earthquakes  in  the  frequency  band  of  roughly  3  to  6  hz,  a 
phenomenon  previously  reported  for  the  same  region  but  with  a  more  limited  data  base. 
This  study  thus  reveals  the  cause  of  the  Pn/Lg  spectral  ratio  discriminant  and  more  impor¬ 
tantly,  the  limitation  of  our  current  theoretical  framework  for  high-frequency  regional 
wave  excitations  for  both  earthquakes  and  underground  nuclear  explosions. 
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1.  INTRODUCTION 

The  seismic  Pn  and  Lg  wave  trains  are  two  of  the  most  prominent  regional  phases  on 
short-period  records  observed  in  continental  areas.  Beyond  a  critical  distance  (generally 
between  100  an  200  km),  the  Pn  and  Lg  waves  become  the  first  and  last  main  high- 
frequency  regional  phases  to  arrive,  untangling  themselves  from  other  prominent  regional 
phases  such  as  Pg  and  Sn.  The  Pn  wave  can  be  modeled  as  an  interference  of  multiple¬ 
diving  waves,  refracted  at  the  Moho  and  traveling  along  the  uppermost  part  of  mantle 
(e.g.,  Cerveny  &  Ravindra,  1971,  Hill,1973).  At  closer  distances,  Pn  is  much  like  a  pure 
head  wave  which  evolves  into  mantle  turning  waves  as  distances  increase  ie.g.,  Sereno  & 
Given,  1990).  At  distances  beyond  about  12-13°,  mantle  turning  waves  are  affected  by 
the  420  and  670  km  discontinuities  and  travel-time  triplication  occurs.  Lg  samples  the 
continental  crust  and  can  be  treated  as  a  sum  of  higher  mode  surface  waves  {e.g., 
Knopoff,  1973),  or  multiple  supercritically  reflected  S  waves  (Bouchon,  1982).  It  can  be 
observed  at  distances  of  up  to  a  few  thousand  kilometers. 

Both  the  Pn  and  Lg  waves  have  been  extensively  used  to  study  earth  structure  and 
seismic  sources.  In  nuclear  explosion  seismology,  they  have  each  been  used  in  magnitude 
determinations  {e.g.,  Nuttli,  1986a,  1986b,  1988;  Vergino  &  Mensing,  1990).  The  ampli¬ 
tude  ratios  of  P/Lg,  including  those  of  Pn/Lg,  have  been  extensively  used  in  studies  of 
discrimination  between  explosions  and  earthquakes  {e.g.,  Blandford  &  Hartenberger, 
1978;  Blandford,  1981;  Nuttli,  1981;  Taylor  et  al,  1988;  Kim  &  Richards,  1993;  Hartse 
etal.,  1997). 

While  the  authors  of  the  studies  of  P/Lg  ratios  generally  agree  that,  at  least  in  certain 
high-frequency  bands,  the  P/Lg  ratios  from  explosions  tend  to  be  higher  than  those  from 
earthquakes,  many  questions  still  exist  on  when  and  why  the  P/Lg  ratio  can  be  used  as  an 
effective  discriminant.  For  example,  it  is  widely  observed  that  the  separation  between 
Pn/Lg  ratios  from  explosions  and  those  from  earthquakes  varies  with  frequency,  paths 
and  source  mechanisms  in  a  complex  manner  {e.g.,  Blandford,  1981;  Chan  et  al.,  1991, 
Lynnes  &  Baumstark,  1991;  Beckers  et  al,  1993;  Kim  et  al.,  1996).  There  appear  to  be 
cases  in  which  Pn/Lg  ratios  observed  for  earthquakes  fall  into  the  explosion  population. 
A  more  profound  problem  is  that,  even  for  the  cases  when  the  Pn/Lg  ratio  in  certain  high- 
frequency  bands  is  known  empirically  to  work  as  an  explosion  discriminant,  there  is  a 
lack  of  understanding  on  why  it  works.  Some  possible  mechanisms  have  been  speculated 
on  based  on  synthetics  {e.g.,  Evemden  et  al.,  1986;  Lilwall,  1988;  Xie  &  Lay,  1994),  but 
it  is  not  practical  to  extensively  verify  these  mechanisms  until  one  can  properly  remove 
the  path  effects  in  the  observed  regional  wave  spectra,  and  obtain  the  spectral  characteris¬ 
tics  of  source  excitation.  Such  removal  of  path  effects  has  been  difficult  and  controversial 
{e.g.,  Mueller  &  Cranswick,  1985;  Harr  et  al,  1986,  Atkinson  et  al,  1997;  Haddon, 
1997)  and  until  recently,  most  studies  on  the  P/Lg  amplitude  ratio  discriminant  have  been 
conducted  by  measuring  the  ratios  themselves;  little  attempt  has  been  made  to  quantita¬ 
tively  infer  spectral  characteristics  of  the  excitation  of  each  phase  (Pn  and  Lg)  by  each 
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source  type  (explosions  and  earthquakes),  to  see  what  might  have  been  the  systematic  dif¬ 
ferences  among  these  excitations  that  cause  the  discriminant  to  work. 

Xie  (1993)  developed  a  non-linear  algorithm  to  simultaneously  invert  for  source 
spectral  parameters,  including  seismic  moment  (Mq),  comer  frequency  (/<-),  and  path  Qq 
and  77  (Q  at  1  hz  and  its  power-law  frequency  dependence,  respectively)  using  seismic  Lg 
waves  from  explosions.  Advantages  of  that  method  include  (a)  it  is  an  exhaustive  search 
method,  thus  requiring  no  starting  model  for  the  unknown  parameters,  (b)  it  allows  Lg  Q 
to  be  path-variable,  and  (c)  it  is  computationally  fast  due  to  a  partitioning  of  model 
parameters,  which  allows  <2o  ^<1  V  ^  be  estimated  with  linear  regressions.  Xie  (1993), 
Xie  et  al.  (1996)  and  Cong  et  al  (1996)  applied  that  method  to  Lg  spectra  from  many 
underground  nuclear  explosions  and  earthquakes  in  central  Asia,  recorded  by  many  mod¬ 
em,  broadband  seismic  stations  in  distance  ranges  of  about  650  km  to  40(X)  km.  They 
obtained  path  Qq,  tj  values  that  correlate  well  with  major  tectonic  features,  and  Lg  source 
Mq  and  /<.  values  for  the  numerous  events.  Based  on  these  values  and  the  respective  val¬ 
ues  for  m^,  the  body  wave  magnitudes,  they  derived  various  scaling  relations  among 
these  values  to  grossly  quantify  the  spectral  characteristics  of  excitation  of  Lg  by  explo¬ 
sions  and  earthquakes.  Among  the  major  discoveries  that  they  made  are:  (1)  for  both 
source  types  the  models  appeared  to  fit  the  Lg  source  spectra  quite  well,  (2)  the  Lg 
Mq  correlates  linearly  with  with  slopes  that  are  either  close  to,  or  slightly  greater  than, 
1.0;  (3)  the  Lg  Mq  tends  to  scale  with  with  a  ranging  from  about  3.6  to  4.0;  and  (4) 
there  appeared  to  be  a  tendency  for  the  explosion  to  be  higher  than  the  earthquake  fc 
at  the  same  Mq  level,  although  there  is  also  a  considerable  overlap  between  the  two  popu¬ 
lations. 

Recently,  Xie  (1998)  presented  a  modification  to  the  method  of  Xie  (1993)  to 
improve  its  stability  in  presence  of  random  and  systematic  errors  in  the  observed  spectra, 
and  successfully  applied  the  modified  method  to  the  1995  western  Texas  earthquake 
sequence  (Af^  -  3.5-5.7).  In  this  report,  we  will  extend  the  method  of  Xie  (1993,  1998) 
to  the  spectral  inversion  of  Pn  waves  to  obtain  estimates  of  path  Pn  Qq  and  rj,  and  Pn 
source  spectral  characteristics.  The  latter  include  the  Mq  and  under  the  idealized  (o~^ 
source  models,  and  the  possible  deviation  of  the  Pn  source  spectra  from  these  models.  We 
will  apply  the  inverse  method  to  Pn  spectra  from  many  seismic  events  in  and  around  the 
Lop  Nor  Test  Site  (LTS)  in  Xinjiang,  China.  These  include  the  last  eight  underground 
explosions  detonated  at  the  test  site  before  September,  1996  (the  time  when  the  United 
Nations  approved  the  Comprehensive  Test  Ban  Treaty),  and  nineteen  earthquakes  that 
occurred  between  1994  and  1996  in  or  near  the  LTS  (Fig.  1).  After  1992,  the  LTS  became 
the  only  active  land-based  test  site  in  the  world,  and  has  provided  all  of  the  new  regional 
seismic  signals  from  underground  nuclear  explosions.  The  large  number  of  stations  and 
high  quality  of  records,  as  well  as  our  new  development  of  the  methodology,  enable  us  to 
derive  grossly  robust  and  reliable  path  corrections  and  source  spectral  parameters  using 
Pn  and  Lg.  In  particular,  for  the  first  time  we  are  able  to  derive  the  scaling  between  Mq 
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and  fc  estimated  using  Pn  from  underground  nuclear  explosions,  and  examine  systematic 
deviations  of  Pn  source  spectra  from  the  idealized  oT^  model.  These  put  tight  constraints 
on  why  the  Pn/Lg  spectral  ratios  may  be  be  used  as  an  explosions  discriminant,  and  why 
the  performance  of  this  discriminant  varies  with  frequenc  y. 


2.  MODELING  OF  Pn  SPECTRA 

Street  et  al  (1975)  were  the  first  to  model  Lg  wave  spectra  via  a  stochastic  approach. 
Sereno  et  al  (1988)  extended  that  modeling  to  the  Pn  wave.  Here,  we  adapt  the  modeling 
of  Sereno  et  al  and  assume  that  A,(/),  the  Pn  wave  spectrum  observed  at  an  ith  station, 
can  be  modeled  as 


A.(/)  =  5(/)/c(^,)G(Ai)exp 


fc(/)r/(/) 


(1) 


where  /  is  frequency,  A„  Vg  and  Qi(f)  are  the  azimuth,  distance,  group  velocity  and 
quality  factor  of  Pn  from  the  source  to  the  ith  station.  /?(0,)  is  the  non-isotropic  source  to 
Pn  radiation  pattern.  X,(/)  and  r,(/)  are,  respectively,  the  site  response  and  the  effects  of 
randomness.  S(/)  is  the  isotropic  component  of  the  the  Pn  source  spectrum,  which  is 
given  by 
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where  Mq  and  are  seismic  moment  and  comer  frequency,  Ps  and  are  source-zone 
density  and  P  wave  velocity,  which  in  this  study  are  set  to  be  2.6  glcm^  and  5.2  km/s  for 
the  explosions  in  the  LTS,  and  2.7  gfcm^  and  6.5  km/s  for  the  earthquakes  (e.g.,  Matzko, 
1992;  Li  et  al,  1995).  The  latter  values  are  also  used  for  the  cmstal  averaged  density  and 
P  wave  velocity,  pc  and  (e.g.,  Roecker  et  al,  1993;  Gao  &  Richards,  1994).  In  Eq.  (2), 
we  have  used  the  omega-square  model  without  overshoot  for  the  earthquakes  (often 
referred  to  as  the  "Bmne’s  model")  and  the  omega-square  model  with  overshoot  for  the 
explosions  (the  modified  Mueller-Murphy  (M.M.M.)  model;  Sereno  et  al,  1988),  with  p 
being  the  parameter  quantifying  the  amount  of  overshoot  (P  is  0.75  for  a  Poisson 
medium).  G(A,)  in  Eq  (1)  is  the  geometrical  spreading  factor  and  takes  the  form 

G(A,.)  =  4S‘(VA,r  (3) 


where  Aq  is  a  reference  distance,  and  m  is  the  decay  rate  of  Aj(/)  at  large  distances 
(A  >  Aq).  Typical  values  for  Aq  and  m  are  about  1  km  and  1.3  to  1.5  for  Pn,  respectively, 
and  are  subjected  to  fairly  large  uncertainties  (Sereno  et  al,  1988;  Sereno  &  Given, 
1990).  In  this  study  we  use  an  m  of  1.3  unless  otherwise  specified.  Qi(f)  in  Eq  (1)  is  the 
apparent  Q  (quality  factor). 
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Qiif)  =  Qoi 


where  j2oi  and  77,  are  Q  to  the  ith  station  at  1  Hz  and  its  power-law  frequency  dependence, 
respectively.  Following  Xie  (1993),  to  develop  a  practical  algorithm  for  the  inversion  of 
Pn  spectra,  we  define  Q'qi  and  rj'i  as  the  apparent  Pn  Q  at  1  hz  and  its  power-law  fre¬ 
quency  dependence,  via 

R(«,>  exp  -  X,(/)  =  exp  -  .  (5) 

V  ygQoi  J  y gUoi  1 


With  this  definition,  Eq  (1)  simplifies  into 


Ai{fj)  =  S(f)0{Ai)exp  -  ‘/7y—  nifj) 


where  we  have  used  subscript  j  to  indicate  the  jth  discretized  frequency.  Equation  (6)  is  a 
simplified,  but  practical  stochastic  modeling  that  relates  the  observed  A,(/,)  to  unknown 
source  and  path  spectral  parameters,  defined  by  a  model  vector 

f  ,  ,  .  .  f 

nJ  ==  Mq, 

where  N  is  the  number  of  stations  recording  the  event,  and  the  size  of  m  is  IN  +  2.  From 
here  on  we  will  drop  the  superscript  in  m  and  it  is  understood  that  estimates  of  0oi,  Voi 
will  be  those  of  apparent  Pn  Q. 


3.  INVERSE  METHOD 


We  now  extend  the  non-linear  method  of  Lg  spectral  inversion  of  Xie  (1993,  1998) 
into  the  Pn  spectral  inversion.  The  inversion  is  event-based,  with  an  objective  to  find  an 
optimal  model,  m,  that  maximizes  <Ta/(ot),  the  posterior  probability  density  function 
(Tarantola,  1987)  which  in  our  particular  case,  may  be  expressed  as 

i=N  j=J{i)  r  -1] 

c'’m(w*)  =  (^onst  X  pjjn)  x  exp  ^  -  Z  2  f  ’ 

where  7(0  is  the  total  number  of  frequencies  available  at  the  ith  station,  r,(/y)  is  defined 
in  Eq  (6)  and  is  calculable  for  any  given  set  of  observed  A,(/y)  and  model  vector,  m. 
Pmi^)  in  Eq  (8)  is  the  marginal  density  function  of  apriori  knowledge  on  the  model. 
Assuming  that  we  possess  apriori  knowledge  on  quantities  (n=l,2,...),  and  that  these 
quantities  are  statistically  independent,  we  have 

n 

where  PxS^)  is  the  apriori  knowledge  on  of  the  form 
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(9) 


PxX”^)  = 


1 


Cn2  ~  C/il 


H(x„  -  xT"  +  c„,)  -  H(x  -  xT^  -  c^)] 


(10) 


with  H{x)  being  a  step  function  and  being  the  apriori  estimate  of  x„.  (x^'^''  -  c„i, 
j^aprtor  ^  defines  an  interval  in  which  x„  may  take  values  with  an  equal  probability. 
Examples  of  c„i  and  c„2  are  (a)  they  equal  each  other  and  approach  zero,  in  which  case 
Pxjjn)  becomes  a  delta  function  and  x„  is  precisely  known,  (b)  c„i  and  c„2  are  the  esti¬ 
mated  error  in  the  estimated  and  (c)  c„i  and  c„2  approach  infinity  and  we  have 

null  apriori  information  on  x„.  Examples  of  x„  include  (i)  /<.,  with  estimated  with 
the  empirical  Green’s  function  approach;  (ii)  <2o^  and  fj,  the  path-average  of  <2o£  77,  : 


Qo=- 


(11) 


ILnAi/Qoi 

i 


The  search  for  optimal  m  that  maximizes  Pmi^)  proceeds  by  a  loop  over  all  possible  Mq, 
fc  values  for  the  event.  In  each  loop,  one  first  estimates  the  path  Qo/,  Vi  values  corre¬ 
sponding  to  the  Afo,  fc  values  using  the  relationship 


1 

— ^  In 

71  Xi 


(  A.(f\ 

^l\J  / 

l,G(A,)5(/) 


J=/^-'''7(2o/  +  e.(/) 


1  =  1,2, 3,... ,7  (13) 


where  e,(/)  is  a  small  residual  to  be  minimized.  Once  Qoi  Vi  are  obtained,  the  com¬ 
plete  set  of  model  vector,  m,  is  checked  to  see  if  it  is  consistent  with  the  apriori  knowl¬ 
edge  in  Eq  (9).  Only  those  m  that  give  non-zero  values  of  Eq  (9)  are  saved.  After  loop¬ 
ing  over  all  of  the  possible  Mq,  fc  values,  the  residual  squares  (the  term  inside  the  expo¬ 
nential  function  in  Eq  (8))  saved  for  all  m  that  are  searched  in  the  loop,  and  give  nonzero 
values  of  Eq  (9),  are  brought  together  and  compared.  The  model  vector  m  that  gives  the 
minimum  residual  squares  (therefore  maximum  Pm(fn))  is  selected  as  the  optimal  solu¬ 
tion.  In  practice,  solving  for  Eq  (13)  requires  a  non-linear  inversion  via  iterative  proce¬ 
dure,  which  slows  down  the  computation.  Therefore,  whenever  the  left-hand  side  of  Eq 
(13)  is  not  negative  for  all  frequencies  at  the  ith  station,  we  take  the  logarithm  of  Eq  (13) 
and  fit  log(0Oi)  (1  ~  Vi)  via  linear  regressions  (c./  Eq  (18)  and  (19)  of  Xie,  1993), 
which  greatly  accelerates  the  computation. 


4.  DATA 


The  Pn  data  set  used  in  this  study  consists  of  vertical  component  Pn  waves  from  the 
last  eight  LTS  underground  nuclear  explosions  and  nineteen  nearby  earthquakes,  recorded 
by  the  broad-band  IRIS  stations  AAK  and  BRVK,  CDSN  stations  (for  earthquakes  only) 
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WMQ  and  LZH,  Kazakhstan  Network  (KZNET;  c.f.  BCim  et  al,  1996)  stations  MAK, 
KURK  and  TLG,  and  the  ten-station  Kyrghizstan  Network  (KNET;  Vernon,  1991).  The 
explosions  occurred  between  September  1992  and  July,  1996,  and  earthquakes  between 
1994  and  1996.  Before  September,  1992  the  station  density  in  the  area  did  not  permit 
multiple-station  recording  of  Pn  at  the  proper  distance  range  (roughly  200  km  ^  A  < 
1300  km).  Lg  spectra  from  many  of  these  events  were  not  studied  before,  and  are  there¬ 
fore  also  inverted  in  this  study  using  the  method  of  Xie  (1993,  1998)  The  sampling  rate 
of  the  instruments  varies  between  about  20  per  second  (the  broadband  channels  of  IRIS 
and  CDSN  stations)  and  100  per  second  (the  high-gain  channels  of  the  KNET  stations 
before  late  1994).  For  explosions  the  CDSN  stations  do  not  provide  data,  making  the 
minimum  distances  of  data  collection  to  be  greater  than  about  790  km.  For  earthquakes, 
data  from  CDSN  stations,  including  station  WMQ,  are  available,  resulting  in  a  minimum 
distance  of  about  200  km.  Between  about  200  and  1300  km,  the  highest  frequencies  at 
which  Pn  spectra  can  be  retrieved  with  sufficient  signal/noise  (S/N)  ratios  (>  2)  are 
dependent  on  the  events’  size,  and  are  typically  between  about  7  and  12  hz.  For  the  Lg 
spectra  recorded  in  the  same  distance  range,  the  respective  frequencies  are  somewhat 
lower  (between  about  5  and  10  hz).  For  Lg  recorded  at  greater  distances  (between  about 
1,400  km  and  2,700  km),  the  highest  frequencies  that  yield  sufficient  (>  2)  S/N  ratios 
range  are  typically  between  about  1  and  5  hz.  For  both  phases  the  lowest  frequency  avail¬ 
able  for  oar  analysis  is  about  0.2  hz. 

A  suiking  feature  of  the  Pn  waves  in  the  study  area  is  the  variability  of  their  ampli¬ 
tudes:  for  the  Lop  Nor  explosions,  Pn  amplitudes  vary  drastically  with  recording  station 
within  the  Kyrghizstan  Network,  a  ten-station  network  with  a  fairly  small  aperture  (less 
than  200  km;  Figure  2).  Figure  3  shows  an  example  of  the  explosion-generated  record 
sections,  where  the  time  domain  peak-to-peak  Pn  amplitudes  vary  by  a  factor  of  19 
between  KNET  stations  KBK  (A  =  1173.7  km)  and  AML  (A  =  1277.5  km).  Figure  4 
shows  the  Fourier  spectra  calculated  using  the  first  4.5  s  of  Pn  in  the  time  series  plotted  in 
Figure  3,  with  a  20%  cosine-taper  window.  It  is  obvious  that  the  amplitude  variations  in 
Pn  spectra  primarily  occur  at  lower  frequencies  (-  1  hz).  For  example,  the  amplitude 
ratio  between  stations  KBK  and  AML  is  about  30  or  greater  between  about  0.8  and  1.0 
hz.  Xie  (1996)  suggested  that  the  most  likely  cause  of  the  amplitude  variation  is  some 
deep-seated  3D  structural  anomaly,  such  as  a  3D  Moho  topography.  In  addition  to  ampli¬ 
tude  variations  of  this  kind,  Pn  amplitudes  from  earthquake  sources  often  show  some 
azimuthal  variations  that  are  not  present  in  Pn  from  explosions.  Specifically,  for  some 
earthquakes  the  ratio  of  Pn  amplitudes  at  the  western  stations  {e.g.,  the  BCNET  stations 
and  station  TLG),  with  respect  to  those  at  stations  in  the  north  (stations  MAK  and 
KURK),  are  systematically  lower  than  the  respective  ratios  for  the  explosions. 

These  variations  make  it  very  difficult  to  conduct  Pn  spectral  inversions  using  the 
stochastic  model  (Eq  (6))  and  the  method  described  in  the  last  section.  To  reduce  the  dif¬ 
ficulty,  for  each  seismogram  containing  Pn,  we  obtained  the  Fourier  spectra  with  a  series 
of  20%  cosine  taper  windows  that  have  a  constant  effective  length  of  4.5  s,  but  increase  in 
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Table  1.  Last  Eight  Lop  Nor  Explosions 


Event 

ID* 

Origin 

Time(H:M:S) 

mb 

Location 
(°N,  °  E) 

No.  of 
Stationsf 

092592 

07:59:59.9 

5.0 

41.763, 88.387 

5 

100593 

01:59:56.6 

5.9 

41.667, 88.695 

6 

061094 

06:25:57.8 

5.7 

41.527,  88.710 

4 

100794 

03:25:58.1 

5.9 

41.662, 88.753 

9 

051595 

04:05:57.8 

5.7 

41.603, 88.820 

12 

081795 

00:59:57.8 

5.5 

41.587,  88.782 

10 

060896 

02:55:57.7 

6.0 

41.657, 88.690 

13 

072996 

01:48:57.8 

4.9 

41.824,  88.420 

10 

The  origin  times,  locations  and  magnitudes  are  from  the 
U.S.  Geological  Survey  preliminary  determination  of  epicen¬ 
ters  (PDE)  bulletin. 

*  Event  ID  is  composed  of  three  groups  of  two-digit  num¬ 
bers,  indicating  the  Month,  Date  and  Year  of  the  event,  respec¬ 
tively. 

t  Number  of  stations  that  provided  Pn  spectra  for  this  study. 
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Figure  1.  Map  showing  locations  of  Lop  Nor  Test  Site  explosions  (stars),  nearby  earthquakes 
(circles)  and  seismic  stations  providing  Pn  and/or  Lg  waveforms  used  in  this  study. 
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Figure  2,  Detailed  map  showing  locations  of  explosions  and  nearby  earthquakes  and  seismic 
stations  providing  Pn  waveforms.  The  great-circle  paths  are  representative  for  Pn/Lg  paths 
from  explosions. 


KNET  Record  Section,  Oct.  7,  1994 
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Figure  3.  Record  section  showing  seismograms  in  the  distance  range  between  1061  and  1298  km,  containing 
Pn  and  Lg  from  the  100794  explosion.  The  first  arrival  is  Pn.  The  straight  line  marks  the  arrival  of  signal  with 
a  typical  Lg  group  velocity  (3.5  km/s).  Numbers  next  to  station  codes  are  peak-to-peak  Pn  and  Lg  amplitudes 
in  digital  counts.  Note  Pn  exhibits  a  large  amplitude  variation. 


MAK,  840.0  km  ULHL,  1060.6  km  KZA,  1149.1  km  KBK,  1173.7  km 


apn|i(duiy 


apn;i|duiv 
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Figure  4.  Fourier  amplitude  spectra  of  the  first  4.5  s  of  Pn  waves  calculated  using  seismograms  from  explo¬ 
sion  100794  (Figure  3).  Note  the  large  amplitude  variations  near  1  Hz  at  varying  stations.  Panel  in  bottom 
right  shows  ratio  of  amplitude  at  KBK  with  respect  to  amplitude  at  AML. 


event  lapse  time  with  a  50%  overlap,  to  successively  cover  Pn  and  early  Pn  coda  with 
group  velocities  (Vg)  greater  than  6.6  km/s.  These  spectra  are  than  averaged  to  reduce  the 
amplitude  variation  associated  with  the  onset  of  Pn.  Averaged  spectra  from  Pn  and  early 
Pn  coda  have  been  previously  used  by  Zhu  et  al.  (1991)  to  stablize  the  Pn  spectral  ampli¬ 
tudes,  and  are  necessary  in  this  study  due  to  the  large  amplitude  variations  in  the  first  few 
seconds  of  Pn.  From  here  on  we  will  refer  to  the  averaged  spectra  of  Pn  and  its  early  coda 
as  "Pn  spectra".  The  Lg  records  that  are  new  to  the  present  study  are  processed  with  the 
previously  described  procedure  (Xie,  1993;  Xie  et  al,  1996). 

5.  SPECTRAL  ANALYSIS  OF  Pn 


§5.7  Spectral  ratios  among  Pn  from  explosions 

The  eight  Lop  Nor  explosions  are  separated  by  a  few  tens  of  km  or  less  (Table  1, 
Figure  2)  and  should  have  very  similar  path  effect  to  permit  an  Empirical  Green’s  func¬ 
tion  analysis  (e.g.,  Li  et  al,  1995).  Here  we  use  the  spectral  ratios,  rather  than  time 
domain  deconvolution,  to  empirically  remove  the  path  effects  and  analyze  the  source 
behavior.  From  Eq  (1),  taking  the  ratio  of  Pn  amplitude  spectra  at  the  same  station  from 
two  colocated  explosion  events  (I  and  II),  we  have 


A'i(f)  S'(f) 
Aj'tJ)  S'Hf) 


with  the  superscripts  denoting  the  event.  The  right  hand  side  of  Eq  (14)  is  the  source 
spectral  ratio  which,  under  the  M.M.M.  model  (Eq  (2)),  can  be  expressed  as 


s'if) 

S"(/) 


M'  [l  +  (1  -  //'  +  //' 

Md  [i + (1  -  2j3)fy  fp  -h  fp 


+  e(f) 


(15) 


where  f  (/)  represents  random,  as  well  as  systematic,  deviations  of  the  spectral  ratio  from 
that  predicted  by  the  M.M.M.  model.  Figure  5  shows  the  station-averaged  spectral  ratios 
among  the  eight  explosions.  There  are  two  interesting  features  in  these  ratios:  First,  the 
events  tend  to  form  two  groups,  each  having  similar  source  spectra.  One  group  of  events 
consists  of  the  two  smaller  (m^  ~  5)  explosions  (092592  and  072996),  and  the  other  one 
consists  of  the  5  to  6  larger  (w^  ~  5.5-6.0)  events.  The  second  feature  in  Figure  5  is  that 
the  ratios  of  five  large  explosions  (M^  ~  5. 5-6.0),  with  respect  to  a  small  explosion 
(072996;  Mf,  =  4.9),  have  a  very  pronounced  overshoot  near  about  2  hz  and  undershoot 
near  4  hz.  Under  the  M.M.M.  model,  Eq  (14)  and  (15)  relate  the  observed  spectral  ratios 
to  the  MqIMq,  //  and  fP  values  of  the  explosions  through  a  non-linear  relation,  which 
allows  us  to  invert  for  these  values  in  a  least-squares  sense  (i.e.,  by  minimizing  the  L2 
norm  of  £(/)).  We  did  such  inversions  by  an  iterative,  linearized  scheme,  with  the  over¬ 
shoot  parameter,  p,  fixed  at  1.0,  a  value  that  is  larger  than  0.75.  The  latter  is  the  expected 
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Figure  5.  Ratios  of  amplitude  spectra  among  Lop  Nop  explosions.  Each  panel  is  for  a 
pair  of  two  events  for  which,  at  each  station,  ratios  of  Pn  spectra  were  calculated. 
These  ratios  are  then  averaged  over  all  available  stations  to  obtain  the  ratios  plotted. 
The  top  two  rows  show  ratios  of  seven  of  the  eight  explosions  studied,  with  respect  to 
the  last  explosion  (072996,  mb  =  4.9).  Smooth  curves  are  theoretical  ratios  based  on 
MMM  model  using  inverted  source  parameters  for  five  explosions  (see  text).  Event  ID 
and  magnitudes  are  listed  for  the  events  whose  spectra  are  used  as  numerator.  Numbers 
inside  the  panels  are  the  best  fit  4-  Pn  from  event  061094  is  recorded  by  only  4  stations 
and  is  excluded  in  the  spectral  ratio  inversions.  The  bottom  row  includes  spectral  ratios 
for  five  larger  events,  with  event  IDs  indicated  on  panel  headers. 


value  for  a  Poisson  medium,  and  is  used  for  the  Lg  spectral  inversion  by  Xie  et  al.  (1996). 
The  MqIMq  ratio  and  //,  //  values  obtained  in  the  inversion  are  used  to  constmct  syn¬ 
thetic  spectral  ratios  (smooth  curves  in  Figure  4),  which  fit  the  overall  trend  of  the 
observed  spectral  ratios,  but  underpredict  the  overshoot  near  2  hz  and  the  undershoot  near 
4  hz.  This  is  because  the  M.M.M.  model,  even  with  the  relatively  large  p  of  1.0,  under¬ 
predicts  the  real  overshoot  in  the  source  spectra  near  the  source  comer  frequencies  (near 
2  hz  for  the  larger  explosions  and  4  hz  for  event  072996,  respectively). 

§  5.2  Estimate  ofPn  Q  using  the  largest  events 

Since  the  Pn  amplitude  spectra  are  unstable  with  varying  paths  and  sources,  it  is 
important  to  have  robust  inversions  such  that  the  resulting  path  Pn  Qq  and  rj  are  similar 
for  the  5  largest  similar  explosions  (event  061094  is  excluded  here  since  only  4  stations 
recorded  it  and  we  are  not  confident  to  infer  if  it  is  similar  to  the  other  larger  explosions). 
We  experimented  with  inverting  Pn  spectra  from  each  event  separately,  with  an  apriori 
knowledge  of  the  form  of  Eq  (9)  and  (10),  such  that  are  only  permitted  to  take  values 
between  1.0  hz  and  3.5  hz.  The  latter  range  is  from  a  conservative  use  of  the  apriori 
knowledge  on  obtained  in  the  last  section  (1.8  to  2.7  hz,  see  Figure  5).  In  the  inver¬ 
sions,  we  obtained  path  Qoi  values  that  vary  significantly.  For  example,  the  Qo  value  esti¬ 
mated  to  station  AAK  is  about  300  when  event  100593  is  used,  and  above  700  when 
event  081795  is  used.  This  indicates  that  the  individual  source  spectrum  deviates  from  the 
M.M.M.  source  iiiociel  in  an  event-variable  way,  and  the  inversion  is  sensitive  to  these 
deviations,  enough  to  cause  the  estimated  model  parameters  to  vary  significantly  among 
the  similar  events.  To  overcome  such  instability  in  model  estimates,  we  averaged  the  Pn 
spectra  among  all,  or  a  subset,  of  the  five  events.  For  each  set  of  event-averaged  spectra, 
we  inverted  for  a  single  m,  which  is  approximately  the  average  model  vector  for  events 
being  averaged.  Owing  to  the  similarity  of  these  events,  this  practice  should  result  in 
robust  and  optimal  estimates  of  Pn  (2o/»  Vi  values.  Figure  6  shows  an  example  of  such  an 
inversion,  in  which  the  observed  Pn  spectra  for  12  broadband  stations  are  averaged  over 
three  events  to  obtain  an  m.  Plotted  in  that  figure  are  the  Pn  source  spectra  calculated 
using  observed  Pn  spectra  with  path  Q  corrections,  and  the  synthetic  source  spectra  con¬ 
structed  using  the  Mq,  obtained  during  the  inversion.  The  overall  fit,  particularly  for 
the  station  average  (last  panel),  is  good  at  high  frequencies,  and  fairly  good  at  lower  fre¬ 
quencies,  except  that  the  synthetic  source  spectrum  has  less  overshoot  than  the  observed. 
This  again  shows  that  the  M.M.M.  model  underestimates  the  level  of  spectral  overshoot 
in  the  real  Pn  source  spectra. 

We  changed  the  event  composition  in  the  averaging  process  by  averaging  two  events 
that  are  not  used  in  Figure  5,  as  well  as  by  averaging  all  five  events.  In  the  respective 
inversions,  the  resulting  path  (2oi  values  to  any  station,  such  as  AAK,  differ  by  less  than 
10%.  The  Mo  values  varied  between  2.6  and  3.2  x  10^^  Nm,  and  fc  values  varied  between 
2.4  and  2.8  hz.  We  then  separately  inverted  each  of  the  five  explosions  again,  by  using  a 
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CHM, 1168, 270  0.6  UCH, 1187, 450  0.2  AAK, 1188, 410  0.4  USP, 1190, 336  0.5  EKS2, 1247, 286  0.5 


Freq.  (Hz)  Freq.  (Hz)  Freq.  (Hz)  Freq.  (Hz)  Freq.  (Hz) 


AML, 1255, 328  0.5  KURK, 1276, 283  0.8  Station  Average 


Freq.  (Hz)  Freq.  (Hz)  Freq.  (Hz) 


Figure  6.  Averaged  Pn  source  spectra  for  three  explosions  (100593,  100794  and 
051595)  at  12  stations.  Fluctuating  curves  are  observed  Pn  source  spectra,  obtained  by 
first  averaging  Pn  spectra  over  the  three  events,  and  then  reducing  these  event-averaged 
spectra  to  source  spectra  by  correcting  for  path  effects.  Path  effects  are  estimated  using 
Qoi,  77i  values  that  are  obtained  by  inverting  event-averaged  spectra.  Smooth  curves  are 
synthetic  source  spectra,  constructed  using  the  MMM  source  model  and  event-averaged 
Mo  and  f,.  values.  The  last  panel  shows  station-averaged  spectra. 
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apriori  knowledge  on  path  Qo/  and  et,  to  individual  stations  defined  by 

=  (//(fir"  +  0. 1 X  fiS”*')  -  //(fir^  -  0- 1  X  2Sf"°")l  (16) 

and 

P,,,(m)  =  [iy(/7r""  +  0.  l)-//(77r''''-0. 1)]  ,  (17) 

where  the  and  values  are  the  average  from  the  three  inversions  of  the  aver¬ 

aged  spectra.  Equations  (17)  and  (18)  permit  the  Qoi  and  rji  to  vary  by  no  more  than  10% 
and  0.1,  respectively,  from  their  apriori  estimates.  These  new  event-based  inversions,  with 
the  apriori  knowledge  on  Qoi  and  77,,  result  in  our  final  estimates  of  Mq  and  fc  values  for 
the  five  largest  explosions. 

As  indicated  in  the  example  in  Figures  3  and  4,  the  estimated  Qq  and  tj  varies  drasti¬ 
cally  across  the  KNET,  due  to  the  large  amplitude  variations  mentioned  before.  The 
^and  ij  (path  averaged  Qq  and  tj  values;  Eq  (11)  and  (12))  are,  however,  quite  robust. 
For  the  eleven  western  stations,  including  the  ten  KNET  stations  and  station  TLG,  Qq  and 
fj  are  364  and  0.5,  respectively.  For  the  two  northern  stations  (MAK  and  KURK)  Qq  and 
fj  are  303  and  0.7.  For  all  paths  the  ^  and  fj  are  354  and  0.5,  respectively,  which  maybe 
taken  as  the  regional  average.  To  see  if  the  inversions  might  have  been  affected  by  the 
abnormal  large  amplitudes  at  stations  KBK  and  KZA,  we  repeated  the  inversions  without 
these  two  stations,  and  found  that  the  resulting  Qq^,  7,  lu  other  stations  are  unchanged. 

Ib  see  whether  the  calculated  Qqi  and  particularly  their  path-averaged  values,  are 
reasonable,  we  compared  them  with  two  estimates  of  path  attenuation  that  are  obtained 
separately.  First,  we  measured  the  two-station  Qq  and  77  between  the  two  northern  sta¬ 
tions  (MAK  and  KURK),  using  available  Pn  spectra  from  two  explosions.  These  two  sta¬ 
tions  are  the  only  ones  that  we  can  find  that  (a)  are  well  separated  in  A;  to  permit  a  rela¬ 
tively  stable  measurement  of  Pn  Q,  (b)  have  similar  azimuths;  and  (c)  are  relatively  free 
of  the  afore-mentioned  drastic  amplitude  variations  among  the  western  stations.  We  cal¬ 
culated  the  two-station  spectral  ratios  from  the  two  available  events,  and  averaged  these 
ratios  to  estimate  the  two-station  Qq  and  77  (c./,  Xie  &  Mitchell,  1990).  Figure  7  shows 
the  result,  where  the  calculated  Qq  and  77  values  are  277  ±  68  and  0.61  ±  0.05,  which 
agree  with,  within  the  estimated  uncertainties,  the  Qq  of  303  and  fj  of  0.7  to  the  two  sta¬ 
tions  (last  paragraph).  A  second  estimate  of  path  attenuation  is  from  Priestley  and  Patton 
(1997),  who  parameterized  the  time-domain  peak-to-peak  Pn  amplitudes  near  station 
WMQ  (Figure  1)  as  decaying  with  distance  at  a  rate  of  ^  distance  of  1,1(X) 

km  and  using  our  modeling  in  Eq  (1)  or  (6),  with  an  m  (Eq  (3))  of  1.3,  that  decay  rate 
would  correspond  to  a  Pn  Qq  of  387,  which  is  close  to  our  estimate  of  Qq  of  354  for  the 
entire  study  area.  We  are  therefore  assured  that  our  (2o,  and  tji  measurements,  together 
with  an  m  of  1.3,  are  grossly  adequate  for  the  paths. 
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.5hz,Qo  =  277  +  -68,  eta  =  0.61  +  -0.05 
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(100794,  051595)  are  recorded  by  the  two  stations  simultaneously,  and  contributed  to  this  plot. 


§  5.3  Sensitivity  of  the  Pn  Q  estimates  to  the  geometrical  spreading  model 

The  estimates  of  Qqi,  77, •  values  in  the  last  section  are,  of  course,  dependent  on  the 
geometrical  spreading  model  assumed.  We  have  chosen  m  to  be  1.3  since  it  is  the  value 
used  in  Scandinavia  by  Sereno  et  al.  (1988),  and  approximates  the  value  of  1.35  at  the 
medium  frequency  used  by  Zhu  et  al  (1991)  for  eastern  Canada.  Xie  (1996)  computed 
extensive  synthetics  using  a  frequency-wave  number  integration  code  and  the  velocity 
models  by  Gao  and  Richards  (1994),  Roecker  et  al.  (1993;  model  Ml)  and  Quin  and 
Thurber  (1992),  with  the  correction  for  the  earth’s  curvature  (Sereno  and  Given,  1990). 
Fitting  the  synthetic  amplitudes  between  0.2  and  2.5  hz  by  a  frequency-independent  geo¬ 
metrical  spreading  of  the  form  of  Eq  (3),  Xie  (1996)  estimated  that  the  m  values  ranged 
between  1.36  and  1.47,  with  considerable  uncertainty  (between  0.2  and  0.3).  To  see  the 
effect  of  varying  m  in  the  spectral  inversion,  we  repeated  the  inversion  using  event- 
averaged  Pn  spectra  with  m  =  1.4.  The  resulting  Mq  and  /<-  values  remained  unchanged, 
but  Qoi,  Vi  values  increased  drastically.  The  fj  for  all  paths,  for  example,  changed 
from  354,  0.5  to  748,  0.3,  respectively.  We  conclude  that  the  Qq  and  77  estimates  in  this 
study  are  valid  only  with  a  frequency-independent  m  of  1.3,  and  may  not  necessarily  rep¬ 
resent  the  true  quality  factor  in  the  uppermost  mantle. 

§  5.4  Spectral  inversions  ofPn  and  Lg  using  apriori  knowledge  on  path  Q 

Generally  speaking,  as  compared  to  the  five  explosions  smaied  in  the  previous  sec¬ 
tions,  Pn  spectra  from  the  remaining  three  explosions  and  ail  of  the  earthquakes  are  more 
affected  by  complications  such  as  fewer  recording  stations  and/or  lower  S/N  ratios  and,  in 
the  case  of  earthquakes,  significant  radiation  patterns.  Spectral  inversion  of  these  events  is 
therefore  conducted  by  using  the  apriori  knowledge  on  path-averaged  Qo  and  v  (Eq  (11) 
and  (12)),  which  is  derived  using  j2of»  Vi  values  obtained  in  the  inversions  of  the  five 
larger  explosions.  To  use  this  knowledge,  we  employed  Eq  (10),  with  the  respective  c„i 
and  c„2  set  at  0  (i.e.,  with  the  right  hand  of  Eq  (10)  taking  the  limiting  form  of  a  delta 
function).  When  Pn  from  the  earthquakes  are  inverted,  the  use  of  the  apriori  knowledge 
on  path-averaged  ^  and  fj  caused  individual  Go/  and  77,  values  to  deviate  significantly 
from  the  respective  values  estimated  in  the  inversion  of  the  five  large  explosions.  For 
example,  the  Go/  values  to  the  northern  stations  MAK  and  KURK,  estimated  during  the 
inversion  of  Pn  from  earthquakes,  are  often  systematically  higher  than  those  obtained  by 
inverting  the  five  large  explosions,  and  just  the  opposite  is  found  to  the  western  stations 
(the  KNET  stations  and  station  TLG).  Lateral  variations  in  Pn  Q  may  have  somewhat 
affected  these  deviations  since  the  earthquake  paths  do  not  exactly  overlap  the  explosion 
paths  (Figure  2).  The  other  likely  cause  for  these  deviations  is  the  effect  of  the  non¬ 
isotropic  radiation  pattern.  Previous  observations  have  suggested  that  Pn  radiation  pattern 
is  a  major  factor  which  affects  the  observed  amplitudes  in  various  areas,  including  Iran, 
Northeastern  U.S.A.,  and  the  study  area  (Nuttli,  1980;  Zhao  and  Ebel,  1991,  Xie,  1996). 
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Using  the  same  methodology  of  Xie  et  al  (1996)  and  Xie  (1998),  we  inverted  the  Lg 
spectra  from  five  LTS  explosions  that  are  not  studied  by  Xie  et  al.  (1996),  and  the  19 
earthquakes.  When  inverting  the  Lg  from  explosions,  the  overshoot  parameter,  fi,  was  set 
at  0.75,  which  was  used  in  Xie  et  al.  (1996)  and  appears  to  fit  the  observed  Lg  spectra 
quite  well.  Figures  8  through  1 1  show  the  resulting  Mo  and  /<-  values  for  both  phases  (Pn 
and  Lg)  and  source  types  (explosions  and  earthquakes),  obtained  in  this  and  previous 
studies  by  Xie  et  al.  (1996)  and  Cong  et  al.  (1996). 

6.  SCALING  BETWEEN  m*  AND  Mq 


Figures  8  and  9  show  the  body  wave  magnitude,  m^,  versus  logio(Mo)  estimated 
using  Pn  spectra  in  this  study,  as  well  as  those  estimated  using  Lg  spectra  in  this  and  pre¬ 
vious  studies.  For  both  source  types  and  both  phases,  the  logarithm  of  Mq  correlates  with 
linearly,  with  slopes  of  about  1.0.  We  conducted  linear  regression  analysis  over  the  m£, 
and  logio(Mo)  values,  and  obtained  the  following  relations; 


log  Mo  =  9. 53(±0. 51)  -I- 1. 16(±0. 09) 
log  Mo  =  9. 26(±0. 45)  +  1. 12(±0. 07) 
log  Mo  =  10. 79(±0. 48)  +  1. 00(±0. 10)  rrif, 
log  A/o  =  9. 96(±0. 48)  +  1 . 17(±0. 08)  ftih 


for  Pn  from  explosions, 
for  Lg  from  explosions', 
for  Pn  from  earthquakes', 
for  Lg  from  earthquakes. 


(18) 

(19) 

(20) 
(21) 


An  interesting  feature  of  these  figures  and  scalings  is  that  for  earthquakes,  at  any  given 
level,  the  Mo  estimated  using  Pn  ("Pn  Mq”)  are  similar  to  Mq  estimated  using  Lg  ("Lg 
Mo").  Both  Pn  Mq  and  Lg  Mq  for  earthquakes  are  higher  than  the  Pn  and  Lg  Mo  esti¬ 
mated  for  explosions  at  the  same  level.  This  is  probably  caused  by  the  fact  that,  at  a 
given  moment  level,  the  amplitudes  of  teleseismic  P  waves  near  1  hz  (the  quantity  used 
for  mi,  determination)  are  systematically  higher  for  explosions  than  for  earthquakes. 


Another  interesting  feature  of  the  figures  and  scalings  is  that,  for  explosions,  at  any 
given  level  the  Lg  Mq  tend  to  be  lower  than  the  Pn  Mq,  although  there  is  some  over¬ 
lap,  making  the  trend  weak.  This  trend  is  probably  related  to  relatively  low  efficiency  of 
Lg  excitation  by  explosions,  as  compared  to  the  efficiency  of  Pn  excitation  by  explosions 
(Sereno  et  al,  1988).  This  trend  should  not  have  been  a  dominant  cause  for  the  Pn/Lg 
spectral  ratio  discriminant  since  it  is  weak,  and  since  it  only  causes  the  ratio  to  vary  at 
low  frequencies  (<  Ihz). 


7.  SCALING  BETWEEN  Mq  AND  f^ 


Figures  10  and  1 1  show  the  versus  Mq  values  obtained  in  this  and  previous  stud¬ 
ies.  Linear  regressions  over  the  logarithm  of  these  values  yield  the  following  relations: 

log  Mo  =  18. 34(±0. 65)  -  4. 73(±1. 2)  log  fc  for  Pn  from  explosions',  (22) 
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explosions.  Straigl 
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earthquakes.  Straight  lines  are  the  linear  regression  fit. 
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this  and  previous  study.  Straight  lines  are  the  linear  regression  fit. 
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this  and  previous  study.  Straight  lines  are  the  linear  regression  fit. 


log  Mq  =  15. 18(±0. 23)  -  3. 81(±0. 39)  log  for  Lg  from  explosions,  (23) 

log  Mo  =  17.08(±0.27)-3.24(±0.28)  log/,  for  Pn  from  earthquakes.  (24) 

log  Mo  =  14. 96(±0. 33)  -  3. 25(d:0. 21)  log  fr  for  Lg  from  earthquakes.  (25) 

For  Pn  from  explosions  (Figure  10;  Eq  (22)),  the  intercept  and  slope  are  subjected  to 
large  uncertainties  due  to  the  small  number  of  points  (8  in  total)  available.  Even  with  the 
uncertainties,  however,  we  may  conclude  that  the  slopes  of  the  above  relations  are  gener¬ 
ally  between  -3  and  -4,  with  those  developed  using  Pn  and  Lg  from  earthquakes  being 
closer  to  -3,  and  those  developed  from  explosions  closer  to  -4. 

The  most  interesting  features  in  Figures  10  and  1 1  are  that  at  the  same  Mq  level,  Pn 
ff.  from  both  explosions  and  earthquakes  tend  to  be  much  higher  (roughly  by  a  factor  of 
5)  than  the  respective  Lg  fr.  For  the  explosions  this  is  an  rather  expected  phenomenon 
since  there  have  been  numerous  studies  reporting  that  the  Pn/Lg  ratios  from  explosions 
are  higher  than  those  from  earthquakes  at  higher  frequencies  {e.g.,  Blandford,  1981;  Tay¬ 
lor  et  ai,  1988;  Kim  &  Richards,  1993;  Hartse  et  ai,  1997).  If  Pn  fr  from  explosions  are, 
but  Pn  fc  from  earthquakes  are  not,  significantly  higher  than  the  corresponding  Lg  fr, 
then  we  may  conclude  that  this  is  the  reason  for  the  Pn/Lg  spectral  ratios  to  work  at  fre¬ 
quencies  higher  than  the  Lg  fr.  What  is  rather  unexpected  is  that  Pn  fr  from  earthquakes 
are  also  significantly  higher  (by  a  factor  of  about  4)  than  the  corresponding  Lg  /,  at  a 
given  Mn  level,  say  Iff^Nm  (Figure  11;  Eq  (24)  and  (25)).  In  terms  of  Mo-  fc  scalings, 
the  discrepancy  between  Pn  /,  and  Lg  fr  from  earthquakes  mimics  the  similar  discrep¬ 
ancy  in  Pn  and  Lg  fr  from  explosions,  and  tends  to  make  both  types  of  sources  similar  in 
terms  of  Mq  -  fr  scaling.  For  earthquakes,  the  amount  that  the  Pn  /,  offsets  the  Lg  fr  at 
the  same  Mq  level  is  about  4,  which  is  slightly  less  than  the  respective  offset  for  explo¬ 
sions  (about  5).  But  this  slight  difference  does  not  seem  to  be  sufficient  to  cause  the 
Pn/Lg  ratio  to  work  as  an  explosion  discriminant  at  higher  frequencies.  It  seems  that  the 
Pn/Lg  spectral  ratio  discriminant  can  not  be  explained  by  a  difference  in  Mq  -  fr  scalings 
for  different  source  types. 

8.  CAUSE  OF  THE  Pn/Lg  RATIO  DISCRIMINANT 

The  scalings  developed  in  the  last  section  enable  us  to  predict,  subject  to  some 
uncertainties,  generic  source  spectra  for  each  types  of  sources  (earthquakes  and  explo¬ 
sions)  and  phase  (Pn  and  Lg),  at  any  given  level.  Figure  12  shows  such  predicted, 
"generic"  source  spectra  for  each  type  of  source  and  phase,  at  the  end  my  levels  of  this 
study  {i.e.,  the  minimum  and  maximum  my  for  each  population  in  Figures  8  and  9).  The 
ratios  between  the  generic  Pn/Lg  source  spectra  are  also  calculated  and  plotted  in  Figure 
12.  If  in  the  observed  data  there  is  a  tendency  for  the  Pn/Lg  spectral  ratio  to  be  higher  for 
explosions  at  certain  frequency  band,  we  ought  to  be  able  to  reproduce  that  tendency,  or 
the  derivation  of  the  scalings  will  have  been  incorrect. 
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Explosion,  mb  =  5.0  Explosion,  nib  =  6.0  Earthquake,  mb  =  4.0  Earthquake, 
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mbS.O,  Eq:mb4.0  Ex:mb5.0,  Eq:mb5.5  Ex:mb6.0,  Eq:nib4.0  Ex:inb6.0,  Eq:mb5.5 


Above  1  hz,  the  Pn/Lg  spectral  ratios  in  Figure  12  all  increase  with  increasing  fre¬ 
quency.  The  manner  of  the  increase,  however,  is  different  for  explosions  and  earth¬ 
quakes.  For  an  explosion,  the  generic  Pn/Lg  ratio  increases  more  rapidly  in  the  frequency 
range  between  the  Pn  /<.  and  Lg  and  reaches  a  maximum  near  the  Pn  due  to  the  Pn 
spectral  overshoot.  In  contrast,  the  generic  Pn/Lg  spectral  ratios  for  earthquakes  increase 
more  slowly  between  Pn  and  Lg  and  no  maximum  is  reached  until  about  10  hz.  As 
a  result,  the  generic  Pn/Lg  ratios  for  explosions  tend  to  be  higher  than  those  for  earth¬ 
quakes  near  the  Pn  This  can  be  more  clearly  seen  in  Figure  13,  where  generic  Pn/Lg 
ratios  for  explosions  are  divided  by  those  for  earthquakes,  resulting  in  "ratios  of  Pn/Lg 
ratios",  which  directly  reflect  the  difference  between  the  generic  explosion  Pn/Lg  spectral 
ratios  and  the  generic  earthquake  Pn/Lg  ratios.  All  of  the  ratios  of  ratios  in  Figure  12 
reach  maximum  near  the  explosion  Pn  /<-  (between  about  2.5  and  5  hz).  This  is  highly 
consistent  with  the  previous  observation  of  Pn/Lg  and  Pg/Lg  ratio  in  the  study  area  by 
Hartse  et  al.  (1997),  who  concluded,  based  on  spectra  collected  from  stations  AAK  and 
WMQ,  that  the  maximum  separation  between  P/Lg  ratios  from  explosions  and  those  from 
earthquakes  occurred  in  a  frequency  band  between  about  3-6  hz.  Thus  the  result  of  this 
study  is  that  the  Pn/Lg  spectral  ratio  from  explosions  are  significantly  higher  than  those 
from  earthquakes  in  a  frequency  band  spanned  by  the  Pn  fc  from  explosions,  mainly  due 
to  the  significant  spectral  overshoot  of  Pn  excited  by  explosions.  We  recall  that  analyses 
in  two  previous  sections  have  indicated  that  the  theoretical  source  model  used  in  this 
study  (the  M.M.M.  model)  underestimates  the  amount  of  spectral  overshoot  in  the 
observed  Pn  from  explosions.  Therefore,  in  Figures  12  and  13  we  have  likely  under¬ 
predicted  the  observed  difference  between  Pn/Lg  spectral  ratios  from  explosions  and 
those  from  earthquakes.  In  other  words,  in  the  observed  Pn,  the  spectral  overshoot  plays 
even  more  significant  role  in  raising  the  Pn/Lg  spectral  ratios  from  explosions  relative  to 
the  Pn/Lg  spectral  ratios  from  earthquakes. 


9.  SOFTWARE  DEVELOPMENT 


The  computer  software  developed  during  the  last  year  includes: 

1.  Program  "nl"  which  takes  observed  data  of  up  to  about  10,000  points,  any  non¬ 
linear  relationship  that  relates  these  data  points  to  up  to  10  model  parameters,  and  finds 
the  model  parameters  that  best  fit  the  data  in  a  least  squares  sense.  The  finding  of  the  best 
model  parameters  is  achieved  by  a  Newton-like,  iterative  scheme.  This  program  can  be 
expanded  into  finding  model  parameters  that  have  larger  size,  with  more  data  points,  with 
minimal  effort. 

2.  Program  "motion"  which  generates  expected  Pn  or  Lg  ground  motion  for  a  given 
type  and  magnitude  of  the  hypothetical  seismic  event  (explosion  or  earthquake),  using 
appropriate  source  spectral  model,  path  Q,  and  source  spectral  scalings  {i.e.,  scalings 
among  m^,,  Mq  and  fc). 
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10.  CONCLUSIONS  AND  DISCUSSION 


Pn  and  Lg  spectra  from  the  last  eight  Lop  Nor  Test  Site  (LTS)  explosions  and  nine¬ 
teen  nearby  earthquakes  are  collected  from  many  broadband  stations  in  the  Republics  of 
Kyrghizstan,  Kazakhstan  and  China  (earthquakes  only).  Pn  from  explosions  show  drastic 
amplitude  variations  (by  a  factor  of  about  20  to  30  in  time  and  frequency  domains)  across 
the  Kyrghizstan  network,  which  is  located  about  1060  to  1280  km  west  of  the  LTS.  Pn 
amplitudes  from  earthquakes  are  further  complicated  by  variations  that  are  not  only  asso¬ 
ciated  with  lateral  variations  of  apparent  Pn  Q,  but  non-isotropic  source  radiation  pat¬ 
terns.  Various  methods  have  been  used  in  this  study  to  improve  the  reliability  in  extract¬ 
ing  the  spectral  characteristics  of  source  excitation  and  path  attenuation  from  the 
observed  Pn  spectra.  We  used  average  spectra  of  Pn  and  its  early  coda  instead  of  the  first 
few  seconds  of  Pn  for  the  analysis.  Averaging  processes  have  also  been  extensively  used 
to  reduce  fluctuations  in  inter-station  and  inter-event  spectral  ratios,  which  lead  to  strong 
constraints  on  Pn  source  spectra  and  path  Q.  The  path-averaged  Pn  Qq  and  77  (Pn  Q  at  1 
hz  and  its  power-law  frequency  dependence,  respectively)  are  obtained  by  inversions  of 
event-averaged  spectra  among  five  large  and  similar  explosions,  and  are  used  as  apriori 
knowledge  for  inversions  of  other  events.  The  improved  methodology,  as  well  as  the 
large  number  of  available  stations,  enable  us  to  invert  for  source  seismic  moments  (Mq) 
and  comer  frequencies  (f^)  using  the  observed  Pn  spectra  and  a  stochastic  model  of  Pn 
excitatioii  and  piopagalion.  Our  anidysis  ucmonstraled  that  Pn  source  spectra  from  explo¬ 
sions  contain  a  spectral  overshoot  that  is  more  significant  than  predicted  by  the  theoreti¬ 
cal,  modified  Mueller-Murphy  source  model,  and  than  the  overshoot  in  explosion-excited 
Lg.  We  also  demonstrated  that  in  the  study  area,  the  estimation  of  Pn  Qq  and  rj  values  are 
highly  dependent  on  the  geometrical  spreading  models  used. 

In  addition  to  analyzing  Pn  spectra,  we  also  inverted  Lg  spectra  from  the  explosions 
and  earthquakes.  Various  scalings  have  been  developed  among  the  body-wave  magnitude 
(m^),  Mq  and  values  estimated  using  both  Pn  and  Lg,  from  both  explosions  and  earth¬ 
quakes.  As  expected,  the  logarithm  of  Mq  values  correlates  with  linearly,  with  slopes 
close  to  1.0.  The  Mq  tend  to  scale  with  /““,  with  a  being  close  to  -3  for  Pn  and  Lg  from 
earthquakes,  and  -4  for  Pn  and  Lg  from  explosions.  One  of  the  most  surprising  discover¬ 
ies  of  this  study  is  that  for  both  explosions  and  earthquakes,  at  a  given  Mq  level,  the  /<- 
values  estimated  using  Pn  are  significantly  higher  than  the  /<.  values  estimated  using  Lg 
for  both  explosions  and  earthquakes.  The  offset  of  Pn  from  Lg  /<-,  at  the  same  Mq 
level,  is  about  a  factor  of  4  and  5  for  earthquakes  and  explosions,  respectively. 

Previously,  it  has  been  reported  that  the  observed  Pn/Lg  ratios  from  explosions  often 
tend  to  be  higher  than  those  from  earthquakes;  the  difference  is  maximized  in  limited 
high  frequency  ranges  (3-6  hz  for  the  study  area).  Through  the  examination  of  the  generic 
spectral  ratios  that  were  computed  using  the  various  scaling  relationships  and  the  theoret¬ 
ical  model  used  in  this  study,  we  found  that  the  dominant  cause  for  the  frequency- variable 
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difference  between  explosion-generated  Pn/Lg  ratio  and  earthquake-generated  Pn/Lg 
ratio  is  the  significant  spectral  overshoot  in  explosion-generated  Pn  spectra.  That  effect 
was  clearly  seen  in  the  generic  Pn/Lg  ratios,  but  should  be  even  more  significant  in  the 
observed  Pn/Lg  ratios  since  the  generic  Pn/Lg  ratios  under-predict  the  effect  of  the  Pn 
spectral  overshoot  for  explosions,  a  limitation  resulting  from  the  modified  Mueller- 
Murphy  source  model. 

The  amount  of  Pn  spectral  overshoot  from  underground  nuclear  explosions  is  likely 
underpredicted  by  currently  available  theoretical  source  spectral  models  other  than  the 
modified  Mueller-Murphy  model.  More  work  is  needed  to  improve  these  models.  The 
theoretical,  stochastic  models  for  Pn  propagation  and  attenuation  are  also  limited  in  that 
they  can  not  precisely  predict  the  effects  of  complicated  velocity  structure.  We  have 
demonstrated  that  even  a  small  amount  of  error  in  the  geometrical  spreading  model  can 
cause  a  large  change  in  the  estimated  Pn  Qq  and  tj  values.  The  estimated  Pn  go  and  t]  val¬ 
ues,  even  when  averaged  over  many  paths,  do  not  necessarily  represent  the  true  quality 
factor  of  the  Earth’s  uppermost  mantle. 

At  individual  stations,  the  observed  Pn  spectra  are  highly  affected  by  the  3D  earth 
velocity  structure.  The  unstable  path  effects,  as  well  as  the  non-isotropic  source  radiation 
patterns  of  earthquake  sources,  make  Pn  spectral  amplitudes  highly  unstable.  Conse¬ 
quently,  the  Pn/Lg  spectral  ratios  are  highly  dependent  on  the  path  and  event  focal  mech¬ 
anism.  The  use  of  average  spectra  of  Pn  and  early  Pn  coda  may  reduce,  but  not  eliminate, 
this  complication,  in  the  worse-case  scenario,  variations  in  Pn/Lg  ratios,  caused  by  vary¬ 
ing  effects  of  path  propagation  and  and  source  radiation  pattern,  may  overwhelm  the  dif¬ 
ference  in  the  Pn/Lg  source  spectral  ratios  caused  by  varying  source  types.  Careful  cali¬ 
bration  of  path  effects  and  increasing  station  density  are  necessary  for  successful  discrim¬ 
ination  of  explosions  from  earthquakes  if  the  spectral  ratios  of  Pn  and  early  Pn  coda,  with 
respect  to  Lg,  are  to  be  used. 

The  findings  of  this  study,  including  the  dominant  cause  of  the  frequency-varying 
performance  of  the  Pn/Lg  spectral  discriminant,  is  limited  by  the  limited  ranges  in  epi- 
central  distance  (>  800  km  for  explosions),  frequency  (usually  between  about  0.2  and  8  to 
10  hz),  and  magnitude  (4.0  to  5.5  for  earthquakes  and  5.0  to  6.0  for  explosions).  It  is  also 
limited  by  the  use  of  only  two  regional  phases  (Pn  and  Lg).  It  will  be  important  to  extent 
the  analysis  to  Pn/Lg  observed  for  other  ranges  of  distance,  frequency  and  magnitude, 
which  may  require  data  sets  that  are  collected  from  a  different  geographic  region.  Other 
important  future  work  will  be  to  extend  the  analysis  to  other  regional  waves,  such  as  Pg, 
whose  ratio  with  Lg  has  been  empirically  suggested  as  being  a  more  robust  discriminant 
between  explosions  and  earthquakes  for  the  study  area  (Hartse  et  al.,  1997). 
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